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Abstract 

The renormalization-group method is used to analyze the low-temperature 
behaviour of a two-dimentional, spin-s quantum Heisenberg ferromagnet. A 
set of recursion equations is derived in an one-loop approximation. The low- 
temperature asymptotics of the correlation length and the uniform suscep- 
tibility are obtained. For small spins (s = 1/2,1) the results are essentialy 

different from those in the spin-wave theory. 
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The revival of the interest in the theory of quantum magnetism has led to creation of 
new approaches of investigation. An effective continuum field theory, which is quantum- 
mechanical generalization of the classical nonlinear cr-model was derived from lattice, large 
spin Heisenberg model of antiferromagnets [|I|]. The two dimensional antiferromagnet on a 
square lattice was treated also by means of the Schwinger boson representation of the spin 
algebra. This representation allows an appropriate mean-field theory of the low-temperature 
behaviour of the sytem ||. At the same time, a modified spin- wave theory of Heisenberg 
(anti)ferromagnets was formulated HJH]. The usual spin- wave theory was suplemented with 
the constraint that the magnetization at each site is zero. This ensures that the sublattice 
rotational symmetry is not broken. 

The results obtained in the papers are in quantitative agreement, but, as was recognized 
by the authors, it is difficult to continue them to smaller values of the spin. 

The aim of the present investigation is to obtain a better understanding of a small 
spin quantum ferromagnet. The renormalization-group method is used to analyze the low- 
temperature behaviour of a two dimensional, spin-s quantum Heisenberg ferromagnet. A 
set of recursion equations is derived in an one-loop approximation. The spin-recursion 
relation shows that the effective spin increases (the spin rescaling factor is two) and hence 
in the limit of infinitely many recursions, the spin-wave approximation can be used. The 
low-temperature asymptotics of the correlation length and the uniform susceptibility are 
obtained solving the recursion relations. For small spins (s = 1/2, 1) the results are essentialy 
different from those in the spin-wave theory. 

The spin-s quantum Heisenberg ferromagnet is defined by the hamiltonian 

h = -J y £§ i - §i (l) 

<i,j> 

where Si are spin operators on site of a two dimentional square lattice with number of 
sites N and lattice spacing a. By < i, j > I denote the sum over the nearest neighbours. 

To preserve the sublattice rotational symmetry in spin-wave theory of the disordered 
phase, Takahashy || imposed a condition of zero sublattice magnetization (Sf) = 0. Using 
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Holstein-Primakoff representation for the spin operators, one can rewrite it as a condition 
that the total number, on average, of spin-waves per site is s. To enfors the constraint, one 
introduces a new term in the hamiltonian h — > h — fJ,J2Sf, with the chemical potential \i to 

i 

be determined from the condition. The system is disordered at any temperature T ^ and 
the chemical potential is positive fx(T) > 0. When temperature goes to zero fx(T) reaches 
and the spin-wave bosons Bose-condense at zero wave-vector. Due to Bose condensation at 
T = 0, the long-range ferromagnetic correlation is present. 

The low-temperature features are determined by the behaviour of the system near the 
origin of the Brillouin zone. Following the renormalization group method one devides the 
Brillouin zone into two parts. Then the half of the zone which contains small wave-vectors is 
of interest to us, and the theory should be reformulated using the reduced Brilloun zone. In 
the site representation, this is equivalent to introduce two sublattices, A and B. I identify 
the sublattices and refer to the new lattice, with number of sites N\ = \N and lattice 
spacing a\ = y/2a , as L\. The hamiltonian (1) can be rewritten in the form 

h = -JEE (t ■ C + t • - V E (Si M + Si**) (2) 

where the sum is over the two directions a x = (a, 0) and a y = (0, a) and over the sites of the 
new lattice L\. The spin operators Si and Si can be realized using canonical bose creation 
and annihilation operators af + , af and af + , af (Holstein-Promakoff representation) Then, 
to leading order of s" 1 the hamiltonian (2) takes the form 

h = sJj2H + ^t+aX+a v + «f + «f + ^a^ +au - af+^a? - a? + a? +av (3) 

-atfaAF - af + af +av ) + n E + *f + «f - **) ~ 

ieLi 

The hamiltonian (3) is diagonalized by the canonical transformation 
d,? = ^(a i + $ i ),a? + = -j=(at + P?), af = ±= - $ t ) , af + = ±= (af - fc) , (4) 

In momentum space the wave-vectors run the half of the Brillouin zone. The dispersion 
of the modes is positive even at T — (fj, — 0) and I refer to them as "fast modes". 

On the other hand, at zero temperature, e^ =Q = 0, and I refer to af and &i as "slow modes" . 
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To proceed further it is more convenient to represent the partition function in terms of 
path integral. In Matsubara representation 



J dfi(a A+ , a A )dfi(a B+ , a B ) exp - j dr £ (a A+ d A + af + af ) + h(r 



(5) 



In formula (5) (3 is the inverse temperature, a i + (t),cl a (t) and a i + (T),af(r) are complex 
functions on imaginery time r, subject to conditions af + (r)af(T ) < 2s and af + {T)af{r) < 
2s, and dfi(a R+ , a R ) is the usual path integral measure for bosons. The fields obey periodic 
boundary conditions. The hamiltonian h(r) is obtained from Eqs (2) replacing the creation 
and annihilation operators by the complex fields. In particular, the spin operators are 
replaced by the vectors Sf(r) and Sf(r) that satisfy S A (r) 2 = s 2 and Sf (r) 2 = s 2 . The 
transformation (4) should be thought of as a change of variables in the path integral (5) 
with Jacobian equal to one. After this transformation, the spin vectors S A and Sj 3 depend 
on the complex fields af(r), a^r) and 0t (r), /3j(r). 

— * 

I shall define the spin-vectors by the equations 



St 



/3+=0,ft=0 



S* 



-5 ? (i)j(o; + , a) 



p+=o,p i= o 

It is convenient to represent the spin-vectors in the form 



(6) 



They depend on the slow modes and satisfy S 2 
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L ft ft q , -A* i g A 

1 5— + fi e< + 



— • 



'1- 



ffff 3 



S (1)i + CpfZi + enpf (7) 



7\ ,„B 



where the real vectors e] = Qj + Qj and = \ (&j — e,) depend on the fields o^~(t) and 

aj(r). They are orthogonal to the vectors S(i)j and to each other and e,- •&j = \ 

The representation (7) is not unique. There is an arbitrariness in the definition of the 

vectors e^, and the coefficients <p A , ip A , <pf, ip B . It is clear from Eqs (7) that a local U(l) 

— * . 

transformation of the vectors and the coefficients does not change the spin vectors Sf and 
Sf. I shall use the following representation for the vectors and 



i 



e 3 = 

^3 



4s 



(4s — ajafj^a 



*t = —(As -at on 



-_ _ J_ + . 
1 4s 



(8) 



The next step of the renormalization group approach is the elimination of the fast modes. 
One can integrate over the fields Pf(r) and Pi{r) only perturbatively. For this reason I 
expand (pf and (pf in powers of (3f{r) and Pi(r) fields. Calculating the coefficients in the 
expansion to leading order of s _1 , one obtains 



rf-y/jft + J^P+ft ^f~_^A+^/3M (9) 

where the terms proportional to Pi Pi are dropped because they do not contribute to the 
effective action. Substituting (9) and (7) in the hamiltonian (2) and selecting the quadratic 
terms in Pf and Pi one gets 

h = ~ E E £ax • h)i + *J E E {pt& + fib + + PtP l ) 

- J7r E E [<%)i • + e^i) A + A + S (1)i ■ (e,a+ + e^) 
-4: E E - ^(Di) 2 (#7% + (io) 



" sJ E E K 1 - + (i - 2e£;) 



f ieLi 

where j = i + a u and the terms proportional to PiPj and Pf Pf are dropped again. 

Integrating over the fields Pf{r) and Pi{r) we lose the rotational invariance of the theory. 
This is well known disadvantage of the spin-wave approximation. To restore the rotational 
symmetry, and to obtain a rotationally invariant effective action one has to average it over 
the elements of the group of rotation || . Under the nonlinear rotational transformations of 
the fields af and oij the scalar product of the vectors (3j and Qj changes the phase e, ■ e 3 - — > 
exp[i(5i — Sj)]Qi ■ Gj, where 5i depends on the fields af and a, and on the group elements. 
Then the module of the scalar product is invariant and the phase transforms as an U(l) 
gauge field on a square lattice 7^ = 7^ + Si — Sj. Hence, the phases enter the invariant 
action in the same way as U(l) gauge fields enter a gauge invariant action. It is not difficult 
to check that these terms have higher dimension and one has to drop them. This means 
to replace <3j • Gj and Gj ■ Gi by |e$ • e_j| in the hamiltonian (10). Straightforward calculation 
using (8) leads to the result 
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3^ (5(i)i 



5, 



(11) 



I shall treat the other terms in the same way. Using the representation (8) and (6), one gets 



S(i)j ■ feat + etOi) + S (1)i ■ fea+ + ejaj) = (S (1)i - S {1)j f + 



(12) 



where the dots stand for terms which become irrelevant after group averaging. 

Integrating over the fast modes I3f(r) and /^(r) only the tadpole graphs contribute to 
the one-loop approximation. Calculating them at zero temperature and taking into account 
the Eqs (11) and (12) I obtain the effective hamiltonian h e s = —J'Y. Yl • 5(i)j> where 
J' = \J 1 — ^ (| + ^2) and j = i + a v . The sum is over the sites of the lattice L\ 
and over the two directions of the previous lattice. I want to represent the sum as a sum 
over the nearest neighbours of the lattice L\. For this reason I write the hamiltonian in 
momentum representation h e g = — J'Yl'S(i)-q • S^i) q (cos aq x + cosaq y ). Let us rotate the 
reduced Brillouin zone q x — (k x + k y ) , q y — ^= (k x — k y ) The new wave- vectors k x and 
k y run over the Brillouin zone of the lattice L\. One can obtain the relation cos aq^+cos aq y = 
\ (2 + cos y/2ak x + cos y/2ak y ^J + 0(k 4 ). Substituting it, and rewriting the expression for the 

— * — * 

effective hamiltonian in site representation one obtains h(i) — — Ji J2 • 

<i,3> 

The effective hamiltonian is a hamiltonian of a Heisenberg ferromagnet with spin 
Si = 2s (see eq (6)), defined on a square lattice with spacing a\ = \/2a, and with exchange 
constant J x = \J [l - £ (§ + ^) . 

Repeating the procedure one obtains the recursion equation 



Sn+i — 2s n , a n +i — V2a n , J n +i — ~~^n 



It is readily to solve them 



2 n s, a n = 22a, J n 



1=1 



1 - 



1 /l 2 

2~T n \2 + ^ 



1 2 

+ - 

it 1 



(13) 



;Jpn{s) 



(14) 



2 l s \2 

I shall consider the equation for the chemical potential. Making the transformation (4) 
and rewriting it in momentum space representation one obtains 
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4- E [<<«,) + OW] = 2s = *i ( 15 ) 

1 g 

where g runs the Brillouin zone of Li lattice. The dispersion of f3 q bosons is positive, 
hence, when the temperature goes to zero, (f3 q f3 q ) goes to zero exponentialy and can be 
dropped in comparisson to (a+a q ). As a result the equation for the chemical potential in 
low-temperature limit takes the form jj-J2{ a t a q) = s i- I n the limit of infinitely many 
recursions, one can calculate < a^a q > in the spin- wave approximation [3]. Then 

-> oo) = ± e -4«8*M8)f> (16) 



where 



p(s) = n 



v " 1/1 2 

~ 2*7 U + 7T* 



(17) 

, =1 L Z" « \Z 7T~ / J 

To obtain the correlation length one has to consider the asymptotic behaviour of the 

two-point static correlation function at long distance. In momentum space this translates 

to small wave-vectors. Each recursion separates smaller and smaller area around the origin 

of the Brillouin zone. In the limit of infinitely many recursions the correlation function 

can be calculated in the spin-wave approximation. Using the relation between the chemical 

potential and the correlation length in spin- wave theory [3], I obtain the low-temperature 

asymptotic of the correlation length 

i 

& - oo) = Jim \ 5 = \afMs^e^^ (18) 

Let us add a new (Zeeman) term in the hamiltonian h — > h — 2HJ2Sf, where H is 

i 

an external constant magnetic field. Then the uniform susceptibility is proportional to 
d 2 J r /dH 2 \ H=ij) where T is the free energy of the system. The recursion does not change the 
magnetic field. For the free energy one obtains T = \T\ + • • ., where T\ is the free energy of 
the system which is defined by the hamiltonian with an additional Zeeman term. The 
dots stand for terms independent on the magnetic field. The multiplier \ in front of T\ is 
a consequence of the fact that the number of sites of the lattice L\ is N\ = \N. Hence 
x(/3 — > oo) = |xi(/3 — > oo) and to leading order of small temperature 
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oo) = Hm — Xn (P 



oo) 



1 




(19) 



37isJp(s) 



where Xn(/3 — > oo) is calculated in the spin- wave approximation [3]. 

From Eq.(17) follows that p(s — > oo) = 1. Substituting this result in formulas (16), (18) 



different from those in the spin-wave theory. 

In conclusion I shall summarize what has been discussed here. Ideas and tools devel- 
oped in the papers of K.G.Wilson §, A.A.Migdal 0, L.P.Kadanoff 0, A.M.Polyakov [§, 
D.Nelson, R.A.Pelcovits, S.Chakravarty and B.I.Halperin |],[K]] have been adapted to for- 
mulate a renormalization group approach in the spin-wave theory. The formalism allows the 
analysis of systems with small spin. 
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and (19), one gets the results from the spin- wave theory [3]. On the other hand p(s)| s= i = 
0.1327 p(s)\ s= i = 0.4462 and one can see that the results for small spins are essentialy 
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